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ABSTRACT 

We study the formation of the spiral structure of barred spiral galaxies, using an 
TV-body model. The evolution of this TV-body model in the adiabatic approximation 
maintains a strong spiral pattern for more than 10 bar rotations. We find that this 
longevity of the spiral arms is mainly due to the phenomenon of stickiness of chaotic 
orbits close to the unstable asymptotic manifolds originated from the main unstable 
periodic orbits, both inside and outside corotation. The stickiness along the manifolds 
corresponding to different energy levels supports parts of the spiral structure. The 
loci of the disc velocity minima (where the particles spend most of their time, in the 
configuration space) reveal the density maxima and therefore the main morphological 
structures of the system. We study the relation of these loci with those of the apoc- 
entres and pericentres at different energy levels. The diffusion of the sticky chaotic 
orbits outwards is slow and depends on the initial conditions and the corresponding 
Jacobi constant. 
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1 INTRODUCTION 

Galactic dynamics together with celestial mechanics have 
played a leading role in the study of orbits as well as in 
the study of order and chaos. Many decades ago and for a 
long time the main interest was about regular orbits since 
they were considered as the main building blocks of galaxies. 
The main idea was that regular orbits are structurally robust 
and therefore they are able to support various morphologi- 
cal features. The most extreme manifestation was elliptical 
gala xies for which no chaotic orbits were con sidered to ex- 
ist (|de Zeeuwlll985l : iBinnev fc Tremalnelll987l ). In rotating 
galaxies the picture was more confused, especially after the 
study of the main periodic orbits in systems representing 
barred galaxies. These studies showed the existence of many 
unstable orbits ( especially near coro tation) which are gen- 
erators of chaos ([Contopoulos 1 1 1 98 lh . Chaos was found in 
the 80 's near the ends of the bar ([Sparke &; Sellwood|[l987l ; 
IContopoulos &; Grosbolll 19891 ). However, chaos was not con- 
sidered as responsible for some of the main features of the 
galaxies like the spiral arms beyond corotation. 

The importance of chaos was realized mainly in the 
90's, with application both in elliptical and in barred spi- 
ral galaxies, mostly in steady state cases (|Merrittl 1996; 
iKaufmann &; Contopoulos! 1 1996h . 

Only in the 2000 's people started to consider asymp- 
totic orbits emanating from the unstable periodic orbits 
in explaining the main features of galaxies, like the spiral 
structure in barred spiral galaxies. The asymptotic orbits 



have initial conditions on the asymptotic manifolds of 
unstable periodic orbits. This study has been done either 
by using analytical potentia l s or by using iV-body sim- 
ulations (IVoglis et al. 
Romero- Gomez et al.l 



2006ai : IContopoulos fe Patsid | 2006| : 
2006, 120071: iTsoutsis et al.1 120081: 



Athanassoula et al.l 



... .200981; lHarsoula Sz Kalapotharakos! 

2009: 1 Athanassoula et aHl2009bh . The former method is not 
self consistent, but it has the advantage of the analytical 
form of the potential, allowing the free choice of all the 
parameters (like, for example, the strength of the bar and 
of the spiral perturbation, or the pattern speed of the bar). 
The latter method has the advantage of self consistency, 
but the time evolving potential and pattern speed imply 
a secular evolution of the system that makes the study 
harder. 

Photometric data have sho wn that spiral structu res are 
related to the old stellar disk ([Grosbol et al.l 12004 ) . Ana- 
lytical models and iV-body dissipationless simulations have 
shown that, the chaotic domains are small in normal spi- 
ral galaxies (non barred), because the spiral perturbations 
are relatively small. Therefore, in many studies up to now it 
was argued that orbits near some stable periodic orbits (reg- 
ula r motions) can support the spiral pattern (see for exam- 
ple iLm fe Shul 11964 lKalnaislll97ll: iLvnden-Bell fe Kalnaid 
Il972l : IContopoulos Ill975l : lToomre!ll977l ). 

In IKaufmann &; Contopoulosl (|l996l ) the authors 
pointed out the role of the so-called 'hot population' 
(|Soarke fe Sellwoodl fl987l ). in supporting the inner parts 
of the spiral arms extending beyond the bar. They found 
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chaotic orbits that wander stochastically, partly inside and 
partly outside corotation. In strong barred galaxies the 
chaotic domains of phas e space have been prov ed very im- 
portant (see for example IShen &; Sellwoodll2004l ). 

The role of the chaotic o rbits in the s piral st ructure 
of a galaxy was pursued by IVoglis et al.1 (|2QQ6bh . using 
self- consistent iV-body simulations of barred spiral galaxies, 
where they found long living spiral arms composed almost 
entirely of chaotic orbits. 

Most chaotic orbits are connected to the asymptotic 
manifo lds o f the main unstable p eriodic orbits. IVoglis et al.l 
(2006a|) and lTsoutsis et all (|2008l ) have considered the apoc- 
entric positions of the asymptotic orbits and the coalescence 
of all the unstable manif olds in a certain range o f energy lev - 
els. On the other hand [Romero- Gomez et al.1 (|2006l . 120071 ) 
and lAthanassoula et al.l If2009al lbh put the emphasis on the 
asymptotic manifolds emanating only from the unstable pe- 
riodic orbits of families originating at the Lagrangian points 
Li and L2. 

In this paper we study a certain snapshot of a 3-D N- 
body simulation of a barred spiral galaxy, in order to give 
a more concrete and explicit description of the mechanism 
of construction of the spiral structure out of sticky chaotic 
orbits which are close to asymptotic orbits from various un- 
stable periodic orbits. Since the iV-body simulation is a self- 
consistent process, it gives all the details of the secular evo- 
lution of the galactic system during a Hubble time, having 
a potential and a pattern speed that vary with time. How- 
ever by using frozen potentials and constant pattern speeds 
that correspond to certain snapshots of the evolution we are 
able to study the role of chaotic orbits in supporting a spiral 
structure that survives for more than 10 bar rotations, i.e. 
for at least 1/3 of the Hubble time. The sequence of these 
snapshots can be considered as an adiabatic approximation 
of the real galaxy evolution. 

Chaotic orbits spend a long time close and along the 
unstable asymptotic manifolds of the unstable periodic or- 
bits in the phase space, due to phenomena of stickiness (see 
IContopoulos fc Harsoulall2008h . In our model, there are en- 
ergy levels where the phase space is not bounded outwards 
and chaotic orbits can escape from the system. However, the 
escape rate is very small, due to stickiness. The orbits spend 
most of their time close to velocity minima on the rotation 
plane of the galaxy. 

We show that the loci of velocity minima are connected 
with the maxima of the density distribution in the configu- 
ration space in every energy level and find the connection of 
these geometrical loci with the apocentres and the pericen- 
tres of the orbits. 

We find the characteristic diagrams of the most impor- 
tant 3-D periodic families and we reveal the spiral structure 
of the galaxy by integrating chaotic orbits that lie close and 
along the main unstable periodic orbits. Finally, we con- 
struct the space configuration by superimposing orbits be- 
longing to different energy levels. The space distribution 
of the apocentres and pericentres of these orbits is very 
similar to the apocentric and pericentric manifolds of 2-D 
asymptotic orbits having initial conditions along the unsta- 
ble asymptotic curves of the 2D periodic orbits. 

The paper is organized as follows: In section 2 we give a 
description of the system providing its main properties. We 
make a frequency analysis separately for the chaotic and reg- 
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Figure 1. The density distribution (in the color scale indicated 
in the top of the panel) of the projected particles of our model on 
the rotation y — z plane made by the superposition of 20 snapshots 
corresponding to the first 20Tf irnct of the integration of the system 
in the fixed effective potential. Black curves denote characteristic 
density contours revealing the spiral structure. Note that the color 
bar shown in this figure applies for all the color scales used in the 
following figures. 

ular populations. In section 3 we discuss the role of the apoc- 
entres and the pericentres of the orbits in comparison with 
the loci of the minimum velocities on the rotation plane. In 
section 4 we discuss the role of the sticky resonant orbits 
in supporting the spiral structure of the system. We present 
examples of specific periodic orbits that are able to recon- 
struct the main features of the galaxy. Moreover, we discuss 
the role of the apocentric and pericentric intersections of 
the 2D asymptotic orbits. In section 5 we give a description 
of the way the orbits escape from the system and finally in 
section 6 we present our conclusions. 



2 DESCRIPTION OF THE MODEL 

The initial c o ndition s of our iV-body model, were created by 
IVoglis et al.1 ([2006b) where four experiments with different 
pattern speeds simulating barred spiral galaxi es were pro- 
duced. In our paper we use the QR4 model of IVoglis et al.l 
(2006b) (the experiment with the greatest value of the bar's 
pattern speed). 

We use 1.5 x 10 5 particles in our simulation. The time 
unit is taken equal to the half mass crossing time of the sys- 
tem (hereafter Thmct)- A Hubble time corresponds to about 
300Thmct- The length unit is taken equal to the half mass 
radius of the system (hereafter r^ m ). Finally, the plane of 
rotation is the y — z plane (intermediate-long axes) and the 
sense of rotation is clockwise. 

The Jacobi constant Ej of an orbit is given by the re- 
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Figure 2. The time evolution of the amplitude of the spiral per- 
turbation 5p/p at r « 1.7rhm- We observe that the spiral struc- 
ture fades out gradually. 
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where V(x, y, z) is the full 3D 'frozen' potential, given by 
a Smooth Field Code (SFC code) as an expansion of a bi- 
orthogonal basis set, x,y,z are the velocities in the rotating 
frame of reference and Q p is the angular velocity of the bar 
(or pattern speed) at the studied snapshot. The value R yz is 

the distance from the rotation axis (^\fy 2 + z 2 ^j . The Jacobi 
constant Ej is called also "energy in the rotating frame" and 
can be written as 



Ej = i(x 2 +y 2 + z 2 ) + V ef f(x,y,; 



(2) 



^QpR yz is the effective 



where V e ff(x,y,z) = V(x,y,z) 
potential. 

The energy of the orbit in the inertial frame is related 
to the Jacobi constant by the following relation 



E — Ej -\- i~2p • L 



(3) 



where L is the angular momentum in the inertial frame of 



reference and fl v 



i being the unit vector along the 



x-axis. Therefore the inertial energy reads 
E — Ej — Q p (yP z — zP y ) 



(4) 



where P y , P z are the momenta in the inertial frame of refer- 
ence. Since 



P y =y + Q P z, 



Vt p y 



(5) 



we have 

E=^(x 2 +y 2 +z 2 )-%(zy-yz) + V(x,y,z) + ^ 2 p R 2 yz (6) 

Our study is made in a distinct snapshot of the iV-body 
simulation that corresponds to 55Thmct- At this snapshot the 
spiral structure of the galaxy is clearly visible and moreover 
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Figure 3. The distribution of the Jacobi constant values for all 
the particles, for the particles moving in regular orbits and for 
the particles moving in chaotic orbits. The maximum of the dis- 
tribution corresponds to chaotic particles around corotation. 



it survives for several rotations of the bar. In order to study 
the mechanism that creates this spiral structure which lasts 
for such a long time we fix the potential and the pattern 
speed and we study the role of chaotic orbits in this fixed 
potential. The corresponding value of the pattern s peed is 
Q p & 27r/8Thmct and is calculated in IVoglis et al.l (|2006b) 
(see Fig. 4a of that paper). In real units this value corre- 
sponds to ~ 25/cra sec -1 kpc -1 . This value is in the range 
of the value s observed in late t ype barred spiral galaxies (see 
for example iFathi et "all l2009h . 

By integrating all the orbits in this fixed effective po- 
tential we reveal the longevity of the spiral structure for at 
least lOOThmct, while after much longer times this structure 
disappears. 

In Fig. 1 we plot, on the rotation plane y — z, the super- 
position of 20 snapshots corresponding to the first 20Thmct 
of the integration time of the system in the fixed effective po- 
tential. Such a superposition consists of more than 2.5 x 10 6 
particles and allows us to reduce significantly the noise that 
appears when the number of particles is small. Since the 
density is now clearly defined, we can safely quantify the spi- 
ral perturbation. In Fig. 2 we show the time evolution of the 
amplitude of the spiral perturbation in our model. More pre- 
cisely, we plot the density excess (Sp/p) caused by the spiral 
structure in a specific radius that corresponds to « 1.7r^ m 
as a function of time. After a time period of lOOThmct, dur- 
ing which the bar has completed ^13 rotations, the spiral 
perturbation is still marginally detectable. This time corre- 
sponds to 1/3 of the Hubble time. 

Figure 3 gives the distribution of the Jacobi constants 
(Eq. [2) for all the particles and separately for chaotic and 
for regular orbits. The distinction between th e two popu- 
lations is done by the method introduced by IVoglis et al.l 
(2006b|). The percentage of chaotic orbits is found to be 
~ 60%. The maximum of the distribution of the chaotic or- 
bits is found between the values Ej(L±) w —1120000 and 
Ej(L±) w —1080000 corresponding to the Jacobi constant 
values of the Lagrangian points and £4,5 respectively. 
From Fig. 3 it is obvious that the regular motions are re- 
stricted inside corotation and therefore they support the 
shape of the bar, while chaotic orbits are responsible for 
the spiral structure outside corotation. 

Figure 4 presents the frequency analysis of the regular 
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Figure 5. The effective potential V e ff along the y-axis (black 
solid line) and along the z-axis (black dotted line). The two curves 
of kinetic energy Ek (dotted and dashed curves) correspond to 
two different energy levels Ej = -1300000 < Ej(Li) (dotted line) 
and Ej = —1050000 > Ej(L\) (dashed line). In the former case 
the minimum velocity corresponds to apocentres for the motion 
inside corotation and to pericentres for the motion lying outside 
corotation. In the latter case the minimum velocity corresponds 
to the position of the maximum of V e f f . 



Figure 4. The frequency analysis of (a) the regular component 
and (b) the chaotic component of the system. The frequency 
ratio q is given by Eq. (7). We see that the vast majority of 
regular orbits lie inside corotation (q > 0) (especially around 
2:1 resonance) while the chaotic orbits lie both inside (q > 0) 
and outside corotation (q < 0). The various spikes shown in (b) 
correspond to chaotic orbits sticking close to various resonances. 



component (Fig. 4a) and of the chaotic component (Fig. 4b) 
of the system. We are only interested in revealing the "disc" 
resonances on the rotation plane and therefore the frequency 
analysis is made in the 2D projection of the orbits on the 
y — z plane. However, there exist also vertical resonances 
that influence mainly the edge-on profiles of our systems 
fsee lPatsis et al.ll2002h . 

The "frequency ratio" q is a parameter given by the 
relation 



^2 ^^T) 



(7) 



where Q is the angular velocity of a particle in the iner- 
tial fra me of reference and k is t he epicyclic (radial) fre- 
quency (Bin nev fe Tremaind Il98/t ). Positive q values cor- 
respond to resonances inside corotation while negative q 
values correspond to resonances outside corotation, which 
are related only to chaotic orbits. In Fig. 4a we see that 
the main resonance of the regular orbits (inside corotation) 
is the 2:1 resonance (or inner Lindblad resonance), while 
there is also a small number of particles around the 3:1 
resonance and an even smaller number around corotation 
(q = 0). The vast majority of these orbits support the bar of 
the system. Figure 4b shows that a rather small percentage 
of chaotic orbits inside corotation lie near important reso- 
nances (e.g. 2:1, 3:1 and 4:1). These represent sticky chaotic 



orbits mainly supporting the out er layers of the bar (see 
lHarsoula &; Kalapotharakosll2009l ) and the innermost parts 
of the spiral structure. The rest of the chaotic orbits move 
around corotation and outside it and some of them stick 
around specific resonances (g=0, -0.5, -1, -1.5, -2, -2.5, -3). 
The q — population corresponds to chaotic orbits located 
around corotation and close t o the PLi, PL2, PL a and PL5 
families (nomenclature after IVoglis et al.ll2006al ). Below we 
will show that the chaotic orbits sticking along the asymp- 
totic manifolds originating from these unstable periodic or- 
bits are responsible for the spiral structure of the system. 



3 THE ROLE OF THE APOCENTRES AND 
PERICENTRES 

In this section we discuss the role of apocentres or pericen- 
tres in revealing the morphological features of the system 
and we relate their geometrical loci with the locus of the 
velocity minima on the rotation plane. 

The distribution of the radial velocities r of the iV-body 
orbits has a preferable concentration around the value zero, 
indicating that the apocentres and pericentres should play 
an important role for the system. 

The geometrical loci of the minima of the velocities v yz 
on the rotation plane should correspond to the maxima of 
the density distribution on the configuration space, because 
particles spend most of their time there. It is of interest to 
investigate how these loci are related to the apocentres or 
the pericentres of the orbits. 

In Fig. 5 we plot the effective potential of a 2D 
approximation of the system (x = 0) along the z-axis 
(Ve//(0, 0, z) = V(0, 0, z) - ^z 2 , black dotted curve) and 
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Figure 6. The distribution of the velocities of the apocentres, of the pericentres and of the velocity minima for (a) Ej < Ej(L±) outside 
corotation (b) Ej < Ej(L±) inside corotation and (c) for Ej > Ej(Li) where the areas inside and outside corotation can communicate. 
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Figure 7. (a) The density distribution of the particles of Fig. 1 when integrated until their (cylindrical) radial velocities become zero 
(R = 0) for the first time, i.e. at the apocentres or the pericentres of their orbits. The density contours corresponding to the initial 
conditions of Fig. 1 are superimposed (black curves), (b) The density distribution of only the apocentres. The comparison with the 
density contours shows that the apocentres support the bar but not the extensive spiral arms, (c) The density distribution of only the 
pericentres. The comparison with the density contours shows that the pericentres support the spiral structure but not the bar. 



along the y-axis (V ef f(0, y, 0) = V(0,y,0) - \Sl%y 2 , black 
solid curve). The maximum of the black dotted curve cor- 
responds to the Jacobi constant Ej(Li) of the Lagrangian 
points Li,2, while the maximum of the black solid curve cor- 
responds to the Jacobi constant Ej^L^) of the Lagrangian 
points Z/4,5- For any Jacobi constant value below Ej(Li), 
for example for Ej — —1300000 (dotted horizontal straight 
line), there is a forbidden region for the orbits (region be- 
tween the two points of intersection with the effective po- 
tential curve). In this case, the phase space, as well as the 
configuration space, consists of two regions that do not com- 
municate with each other (the one inside and the other out- 
side corotation). By plotting the corresponding curve of the 
kinetic energy (dotted curve) and because of (Eq.[2]), we see 
that for the region inside (outside) corotation the velocity 
minima correspond to the apocentres (pericentres) of the or- 



bits. However, for Ej values above Ej(L^)^ for example for 
Ej — —1050000 (dashed horizontal straight line), the region 
inside corotation can communicate with the region outside 
corotation and the minimum velocity Vmin (minimum of the 
dashed curve) is found at the position of the maximum of the 
effective potential (i.e. the minimum difference of Ej — V e ff). 
In this case, the geometrical locus of the velocity minima on 
the configuration space is not directly related to the apoc- 
entres or the pericentres of the orbits. 

Figure 6 presents the distribution of the velocities of 
the pericentres of the apocentres and of the velocity minima 
for three cases: (a) for orbits trapped outside corotation 
with Ej < Ej(Li), (Fig. 6a), (b) for orbits inside corotation 
with Ej < Ej(Li), (Fig. 6b) and (c) for orbits that can 
explore the whole permissible phase space inside and outside 
corotation with Ej > Ej(Li), (Fig. 6c). 
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Figure 8. The density distribution of the particles of Fig. 1 (with 
Ej < Ej{L/Cji) at the positions corresponding to their local plane 
velocity minima (v yz = v m i n ). We observe that this distribution 
fits well the morphology of the galaxy (black contour lines). In 
this figure we have omitted the particles with Ej > Ej{L^) which 
are gathered along the geometrical locus of the maximum of the 
effective potential (gray elliptical curve). 

In Fig. 6a we observe that the pericentres of the orbits 
outside corotation have smaller velocities than the apocen- 
tres. Moreover the distribution of the velocity minima gives 
good agreement with the distribution of the velocities of the 
pericentres. At the same time, the distribution of the veloc- 
ities of the apocentres almost coincide with the distribution 
of the velocity maxima (not plotted in the panel). Therefore, 
the particles outside corotation spend most of their time 
near the pericentres of their orbits which is consistent with 
the description of Fig. 5. On the other hand the opposite is 
true for particles inside corotation (Fig. 6b). The particles 
inside corotation have smaller velocities at their apocentres 
than at their pericentres. Moreover, the apocentric veloc- 
ity distribution almost coincides with the distribution of the 
velocity minima. Here again the distribution of the velocity 
maxima (not plotted in the panel) almost coincides with the 
distribution of the velocities of the pericentres. 

For the orbits having Jacobi constant above Ej(L\) 
(Fig. 6c) even though there is not a clear breakoff of the 
distributions, the pericentres have statistically smaller ve- 
locities than the apocentres. Both pericentres and apocen- 
tres present a peak around the same value of velocity, which 
also coincides with the peak of the velocity minima distri- 
bution. The distribution of the apocentres, however, has a 
second important peak around a somewhat greater value of 
velocity. 

In Fig. 7a we plot, in color scale, the density distribution 
of the particles of Fig. 1 taking each of them at a position 
corresponding to r = 0, i.e. at the apocentre or the pericen- 
tre of its orbit. The black curves correspond to the density 
contours of the space distribution shown in Fig. 1. Figure 7b 



presents, the density distribution of only the apocentres of 
the orbits. The corresponding density maxima inside coro- 
tation support the bar (since they correspond to loci of min- 
imum velocities, see Fig. 6b) as well as a small extension of 
the bar to the inner parts of the spiral structure. However, 
the density maxima outside corotation, do not correspond 
to the density contours of the real spiral arms (black curves) 
but they extend far beyond. Figure 7c presents the density 
distribution of only the pericentres of the orbits. We see that 
these density maxima support a spiral structure very close 
to the real one, but a little more tight, and they do not 
support the bar. 

Thus, we conclude that, in general, the apocentres of 
the orbits support the shape of the bar and the origins of 
the spirals near corotation, while the pericentres of the or- 
bits better support the spiral structure near and outside 
corotation. 

Figure 8 presents, in gray scale, the density distribu- 
tion of the particles of Fig. 1 at positions corresponding 
to local minimum values of their rotation plane velocities 
(v yz = Vmin)- The comparison with the density contours 
(black curves) of the initial conditions shows that the loci of 
the velocity minima reveal satisfactorily the features of the 
galaxy, i.e. the bar and the spiral arms. In this figure, how- 
ever we have removed the overdensity of the particles along 
the theoretical locus of Vmin (gray elliptical curve) that rep- 
resents the maximum of the effective potential (V e ff) for 
Jacobi constants above Ej{L^) : because it does not corre- 
spond to the distribution of the real particles. 

Here we note that a density maximum in some area can 
be produced either by frequent (but not too short-lasting) 
transits or by long-lasting (but not too rare) transits of 
the particles through this area. For Jacobi constant values 
higher than Ej(L±) the locus of the local velocity minima 
(red curve in Fig. 8) corresponds to velocity values that can 
be comparable to the mean velocity of the whole orbits. In 
Fig. 9 we have plotted the density distribution of the par- 
ticles with Jacobi constant values Ej > —977000 > Ej(L±) 
corresponding to low and high values of the ratio A =< 
Vmin > / < Vyz > as indicated on the two panels of Fig. 9. 
We observe that the particles with low A values (Fig. 9a) 
form a ring-like density maximum area lying close to the 
locus of the maximum V e ff (red curve). However the dis- 
tribution is not uniform along the red curve, but it follows 
the distribution of the inner parts of the spiral arms. The 
maximum density area of the particles with high A values 
(Fig. 9b) is located inside the red curve corresponding to the 
velocity minima. In this case the density maxima are due to 
the more frequent transits of the orbits from this area. 



4 STICKY CHAOTIC ORBITS 

In Figure 4 the frequency ratio q of the chaotic population 
is plotted showing preferable concentrations around specific 
resonant periodic orbits. In order to study the behavior of 
sticky resonant orbits we find the characteristics of the most 
important periodic orbits in the 3D case. Figure 10 presents 
the characteristics of the main 3D periodic orbits, i.e. the y 
(z for the PL4 family) component as a function of the Jacobi 
constant. We must point out that the periodic orbits PL\ 
and PL2 exist only above the value Ej(Li) = —1120000, 
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Figure 9. (a) The density distribution of the particles having Ej > —977000 > Ej{L 4 ) and ratio A 



< 0.25 together with 



<Vy Z > 

the locus corresponding to the effective potential maxima V(eff) m ax (red curve). The density maxima in this case are located near the 
curve of V(ef f)max but along particular directions, (b) The same as in (a) but for particles having A = <^ min > > 0.6. The density 

<V yz > 

maxima are located inside the curve of V(eff) m ax ( see text for details). 
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Figure 10. (a) The characteristics of the main 3D periodic orbits of the system, (b) A focusing in the inserted frame of (a). 



while the periodic orbits PL4 and PL5 exist above the value 
Ej(L 4 ) = -1080000. We will show below that all the im- 
portant resonances, presented in the characteristic diagram 
(Fig. 10) and in the frequency analysis of the chaotic orbits 
(Fig. 4b), contribute to the formation of all the morpholog- 
ical features appearing in each Jacobi constant value. 

Figure 11 presents, in color scale, the density distri- 
bution of particles on the rotation plane for nine different 
Jacobi constant levels. The top left panel corresponds to 
the lowest values of Jacobi constants Ej (close to the po- 
tential well) while the bottom right panel corresponds to 



the highest values of Ej*8. The corresponding values of the 
Jacobi constant are marked on the top of each panel. Note 
that outside corotation (in energy levels with Ej > Ej(Li)) 
there exist practically only chaotic orbits (see Fig. 3). A gen- 
eral trend is that the space distribution of particles having 
lower values of Jacobi constants supports mostly the outer 
parts of the spiral arms as well as the shape of the bar, while 
particles having greater values of Jacobi constants support 
the inner parts of the spiral arms. This is obvious in Figs. 
11 b-h, where we see a gradual shift of the maximum of the 
density distribution of particles outside corotation from the 
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Figure 11. The density distribution on the y — z plane of the particles belonging to nine different Jacobi constant bins. The corresponding 
values are marked on the top of each panel. The top left panel corresponds to the smallest values of Jacobi constants Ej (close to the 
potential well) while the bottom right panel corresponds to the highest values of Ej's. In this figure we reveal which structures are 
supported by the particles corresponding to the different Jacobi constant values. 

(|Contopoulos &; Harsoulall 20081 ) of chaotic orbits around the 
unstable asymptotic manifolds of the unstable periodic or- 
bits, in the phase space. We show that this stickiness sup- 
ports the formation of the spiral structure, while the (slow) 
diffusion along the paths of the manifolds is responsible for 
the fadeout of this structure. 

In the next subsections we give four examples of 3D 
chaotic orbits lying near and along some main unstable pe- 
riodic orbits, having initial conditions inside corotation, near 
corotation and outside corotation. 



outer parts of the spiral arms to the inner parts and the 
innermost extensions of the bar (Figs. llf,g). On the other 
hand there are levels of the Jacobi constant that do not sup- 
port at all the spiral structure. This happens for the lowest 
values of the Jacobi constant near the potential well (Fig. 
11a), where the distribution of particles outside corotation is 
quite uniform, and for the highest values of the Jacobi con- 
stant (Fig. Hi), where only chaotic orbits exist, which do 
not support the spiral structure but rather the outer parts 
of the bar. 

Below we give some examples of the effect of stickiness 
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Figure 12. (a) The density distribution on the y — z plane of 20000 3D orbits integrated for several T^mct having initial conditions close 
to and along the unstable periodic orbit —1:1 and Jacobi constant Ej = —1252560. Black curves indicate the high density contours 
of real particles corresponding to the same energy level. We note the coincidence of these density contours with the density maxima of 
the integrated orbits, (b) the density distribution of the pericentres superimposed with the high density contours of the corresponding 
energy level, (c) the same as in (b) but for the apocentres. We observe that the pericentres (apocentres) fit (do not fit) the density 
maxima of the real particles. 



4.1 The -1:1 family 

Figure 12 gives an example of the density distribution of or- 
bits around an unstable periodic orbit of the —1:1 family 
with Ej — —1252560 (within the range of the top-right panel 
of Fig. 11) where the motion is limited outside corotation. 
We have taken 20000 initial conditions all along this 3D un- 
stable periodic orbit (and its symmetric with respect to the 
origin), with small deviations from it, on the 6-dimensional 
space (x,y,z,v x ,v y ,v z ). Then we integrated all these orbits 
for a time corresponding to many rotations of the bar. In 
Fig. 12a we plot the density distribution of these orbits on 
the rotation plane. In this figure we have considered a time 
interval within 5-6 bar rotations. In the same figure we have 
superimposed the high density contours of the correspond- 
ing Jacobi constant (black curves) and we see that the orbits 
close to the —1:1 periodic orbit reproduce well the density 
distribution giving the outer parts of the spiral structure. In 
Fig. 12b (12c) we plot the density distribution of the pericen- 
tres (apocentres) of the orbits of Fig. 12a together with the 
high density contours (black curves) corresponding to this 
Ej value. On the one hand the pericentric density maxima 
seem to follow a part of the spiral structure but at a slightly 
smaller radius. On the other hand the apocentric distribu- 
tion does not fit well the real density maxima. Finally we 
remark that the loci of the velocity minima practically co- 
incide with the ones of the pericentres. 



4.2 The 3:1 family 

Figure 13 gives a similar information as Fig. 12, but for an 
unstable periodic orbit of the 3:1 family, which is located in- 
side corotation for Ej = —1090000 (within the range of the 



middle-right panel of Fig. 11). This Jacobi constant value, is 
close to Ej(Li), and the areas inside and outside corotation 
can communicate. We considered again 20000 initial condi- 
tions all along this 3D (3 : 1) unstable periodic orbit (and its 



symmetric with respect to the origin) , with small deviations 
from it, in the 6-dimensional space (x,y,z,v x ,v yi v z ). These 
orbits follow, in the phase space, the unstable directions of 
the asymptotic manifolds originating from the periodic or- 
bit. For a short time (the shortness depends on how large is 
the initial divergence from the periodic orbit) these orbits 
stay very close to the periodic orbit, but later they start 
deviating considerably from it. In Fig. 13a we see the den- 
sity distribution (in color scale) of these orbits after within 
a time interval of 3-4 bar rotations. The black curves indi- 
cate the high density contours of the real particles for the 
same Jacobi constant value. It is evident that these two dis- 
tributions fit o ne ano t her ve ry well. A similar behavior has 
been found bv lPatsid |2006) for chaotic orbits near the 4:1 
resonance. This type of orbits was considered responsible 
for the inner parts of the spiral arms in an analytical model 
representing the spiral galaxy NGC 4314. 



In Fig. 13b we have plotted the positions of all the 
orbits for the time interval of Fig. 13a. The color of each 
point shows the corresponding velocity value on the rota- 
tion plane. The color scale is the same as for the densities 
(see Fig. 1) which means that blue color corresponds to the 
minima and red color to the maxima of the velocity values. 
We observe that the density maxima are well correlated with 
the areas of small velocity values (dark to light blue colors 
in Fig. 13b). Figures 13c, d are similar to Figs. 12b, c. We 
observe that the pericentres of the orbits (Fig. 13c) form a 
slightly more closed spiral structure than that of the real 
particles (black density contours). On the other hand the 
apocentres of the orbits (Fig. 13d) fit well only the density 
contours of the bar and the innermost parts of the spiral 
structure that originate from the end of the bar. In Fig. 13e 
we plot the density distribution of the loci of the local veloc- 
ity minima. We notice that this distribution matches better 
the distribution of the apocentres in what concerns the bar 
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and the distribution of the pericentres in what concerns the 
spiral structure. 



4.3 The 2:1 (xi) family 

The 2:1 (or x\ after t he nomenclature of 
IContopoulos &; Papavannopoulosl [l980) family is gen- 
erally considered responsible for the formation and the 
robustness of the bar. This is due to its morphology 
together with the fact that this family is usually stable for 
low values of the Jacobi constant. According to this point of 
view the bar should end not beyond the point where the x\ 
family turns from stable to unstable. Here we point out the 
usefulness of the unstable part of this family which supports 
the spiral structure. Figure 14a, b is similar to Fig. 13a, b 
but for the x\ family and for the same Jacobi constant 
value. We observe that the density of the orbits starting 
close to the x\ periodic orbit is similar to the density of 
the orbits with initial conditions near the 3:1 periodic 
orbit. Their density maxima support the morphological 
structures (bar and spiral) of the real particles for the same 
Jacobi constant value (Fig. 14a). These maxima lie near 
the areas corresponding to low values of plane velocities v yz 
(Fig. 14b). 



4.4 The PLi, 2 families 

Figure 15 is similar to Fig. 14 but for the unstable periodic 
orbits PL 1,2 which are located near corotation and for the 
same Jacobi constant value. In this case we establish the 
same behavior (as for the cases 3:1 and x\) for the orbits 
initiating close to the symmetric periodic orbits PL\^- The 
orbits starting near the periodic orbits 3:1, 2:1, PL 1,2 follow 
the unstable directions of the asymptotic manifolds originat- 
ing from these periodic orbits. For the same Jacobi constant 
value these manifolds cannot intersect each other, and this 
means that these manifolds are parallel in the phase space. 
Thus, the orbits along the unstable directions of the various 
manifolds follow parallel paths and present a similar behav- 
ior in the configuration space. The major difference between 
the different sets of orbits starting near the different unsta- 
ble periodic orbits (3:1, 2:1, PLi^) is the diffusion rate (see 
section 5), which increases towards higher resonance orders 
(2:1^3:l^oo:l = PLi, 2 ). 



4.5 Other families 

All the other families of periodic orbits have similar behav- 
ior to those studied in Figs. 12-15. Thus, the behavior of the 
orbits with —1300000 < Ej < El 1 lying outside corotation 
(e.g. —2 : 1) resembles that of the —1:1 family (Fig. 12) 
supporting the outer parts of the spiral structure. Families 
lying inside corotation with Ej > El 1 (e.g. 4:1, PL^s) be- 
have similarly to the families 2:1, 3:1, PLi,2 (Figs. 13-15) 
and support the bar and the inner parts of the spiral arms 
connected to the bar. Note that there are chaotic orbits as- 
sociated with the 2:1, 3:1 and 4:1 families, located inside 
coro tation for Ej < Eli that support the envelope of the 
bar (Har soula &; Kalapotharakos] l2009h . Some of the char- 
acteristics of these families are plotted in Fig. 10. Moreover, 



in Fig. 4b we see many real particles that stick around var- 
ious resonances both inside and outside corotation (see the 
various spikes corresponding to different resonances). 

4.6 Reconstruction of the galaxy by the 
superposition of the main families 

The chaotic orbits close to the unstable periodic orbits corre- 
sponding to different Jacobi constant values may form dif- 
ferent density distributions. However, the orbits at every 
Jacobi constant support the density distribution of the real 
particles at the same value of the Jacobi constant. Here we 
demonstrate what is the contribution of a family, all along 
its characteristic, to the global morphology of the galaxy. 
This is done in Fig. 16 for the 3:1 family. In Fig. 16a we plot 
100 periodic orbits (black points) of the 3:1 family equally 
sampled along its whole characteristic together with the den- 
sity contours of the galaxy (black lines). These orbits are 
located inside corotation. In Fig. 16b we plot the density 
distribution of 20000 orbits having initial conditions close 
and along the periodic orbits of Fig. 16a integrated for 3-4 
bar rotations. Although the initial conditions are all inside 
corotation supporting the bar shape, after a short transient 
period of stickiness along the periodic orbits they are dif- 
fused outside corotation and form a spiral structure which 
is slightly more closed than that of the galaxy (Fig. 16b). 

Figure 17 is similar to Fig. 16 but for the 2:1 (or xl) 
family. Note that the superposition of the periodic orbits do 
not cover the inner parts of the bar since we have chosen only 
the unstable periodic orbits along the characteristic of the 
2:1 family. In Fig. 17b we see that the orbits starting close to 
the periodic orbits of Fig. 17a form a spiral structure fitting 
well the inner part of the spiral structure of the galaxy (black 
curves) . 

Finally, Fig. 18 provides the same information as 
Figs. 16,17 but for the —1:1 family. It is evident that the 
density maxima of the density distribution of these orbits 
support the outer parts of the spiral structure of the galaxy. 
The same picture is also true for the other families outside 
corotation that have q < (e.g. the -2:1 family). 

In Fig. 19 we plot the superimposed density distribu- 
tions of the orbits shown in Figs. 16b, 17b, 18b (color scale) 
together with the high density contours of the galaxy (black 
curves). We observe that the galaxy is well reconstructed 
by this set of orbits. This figure clearly underlines the im- 
portance of the unstable periodic orbits and of the chaotic 
orbits associated with them in forming the spiral structure. 

4.7 2D Asymptotic Orbits 

The main results of the above study, for the 3D orbits along 
specific important resonances of the system can be also ac- 
quired if we study the space distribution of the pericentres or 
apocentres of the asymptotic orbits in a 2D approximation 
of the system (x = 0, x = in Eq.[2]). Namely, if we consider 
2D unstable periodic orbits (the most important for every 
energy level) and take initial conditions along their unsta- 
ble asymptotic curves, the integration of these asymptotic 
orbits can also reveal the morphological features of every 
energy level. In Fig. 20a, 20b we give an example of the 
projection on the configuration space of the apocentric and 
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Figure 13. (a) The density distribution on the y — z plane of 20000 3D orbits having initial conditions close to and along the unstable 
periodic orbit 3:1 (and its symmetric with respect to the origin) at Ej = —1090000. The plotted distribution corresponds to a time 
interval within 3-4 bar rotations. Black curves denote density contours of the real particles corresponding to the same Jacobi constant 
value, (b) The positions of all the orbits at the snapshot shown in (a). The color of each point indicates the corresponding velocity 
value on the rotation plane. The color scale is the same as for the densities (see Fig. 1). (c) The density distribution of the pericentres 
superimposed with the density contours of the corresponding Jacobi constant value, (d) The same as in (c) but for the apocentres. (e) 
The same as in (c) but for the local velocity minima v m i n . 
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Figure 14. Similar to Figs. 13a, b but for the unstable periodic orbit 2:1 (x±) and for the same Jacobi constant value. 
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Figure 15. Similar to Figs. 13a,b but for the unstable periodic orbits PL\^ and for the same Jacobi constant value. 
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Figure 16. (a) The projection on the y — z plane of 100 periodic orbits (black points) of the 3:1 family equally sampled along its whole 
characteristic together with the density contours of the galaxy (black lines), (b) the density distribution on the y — z plane of 20000 
orbits with initial conditions close and along the periodic orbits shown in (a) after integration for 3-4 bar rotations. 
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Figure 17. Similar to Fig. 16 but for the unstable part of the 2:1 family. 
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Figure 18. Similar to Fig. 16 but for the —1:1 family. 
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Figure 19. The superposition of all the distributions shown in 
Figs. 16b, 17b, 18b superimposed with the density contours of the 
galaxy (black curves). All the main morphological features of the 
galaxy are reconstructed by the superposition of these orbits. 



pericentric intersections (black points), respectively, of the 
2D manifold of the —2:1 unstable periodic orbit corre- 
sponding to Ej = —1135000 < Ej(Li) where the areas in- 
side and outside corotation cannot communicate. Note that 
the phase space in the area outside corotation, in the 2D ap- 
proximation of the system, presents no tori or cantori with 
small holes, but only a chaotic sea with escapes to infinity. 
In each panel of Fig. 20 we have superimposed the density 
distribution (in color scale) of the real particles at the same 
Jacobi constant value. We observe that the apocentres do 
not coincide with the maxima of the density while the peri- 
centres trace well the areas of the density maxima, but for 
slightly smaller radii. This is compatible to the information 



provided by Fig. 12 according to which for Jacobi constant 
values below Ej(L±) and outside corotation the pericentres 
reveal the density maxima since they lie near low values of 
the plane velocity. 

Figure 21 is similar to Fig. 20 but for a different Ja- 
cobi constant i£/=- 1120000, which is just above Ej(L±), and 
therefore the area inside corotation has just communicated 
with the area outside corotation. In this energy level and in 
the 2D approximation there exist unstable periodic orbits 
of the PLi, PL2 and —1:1 families. In Fig. 21a,b we plot 
the apocentres and the pericentres, respectively, of all the 
asymptotic orbits (having initial conditions on the unsta- 
ble asymptotic curves) of all the periodic orbits mentioned 
above, superimposed on the density distribution of the real 
particles in the specific Jacobi constant value. In this case 
we see that the pericentres together with some parts of the 
apocentres (mostly of the PLi, PL2 families) trace well the 
density maxima of the density distribution. This is com- 
patible to the information provided by Fig. 13 according to 
which the pericentres and the apocentres for Jacobi constant 
values above Ej(Li) reveal the density maxima when they 
lie near low values of the plane velocity. 



5 RAYS, BELL-TYPE CURVES AND 
ESCAPING ORBITS 

The stickiness effect along the asymptotic manifolds of the 
unstable periodic orbits can last for a number of dynam- 
ical times reinforcing the spiral structure, as it has been 
described in the previous sections, but then the chaotic or- 
bits can escape to very large distances and finally they may 
escape to infinity. This diffusion is relatively slow and its 
time scale can be compared to the age of the Universe. In 
this section we study the way these chaotic orbits escape to 
infinity. 

Figure 22 presents the unstable asymptotic curve of the 
PL2 orbit on a Poincare surface of section z, z for y = 
and y > of the 2D approximation of the system and a Ja- 
cobi constant close to Ej(Li). We have taken a small initial 
segment of 50000 points and of length 10 -6 from the un- 
stable periodic orbit and along the direction of the unstable 
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Figure 20. (a), (b) The projection on the configuration space of the apocentric and pericentric intersections (black points), respectively, 
of the 2D manifold of the —2 : 1 unstable periodic orbit corresponding to Ej = —1135000 < Ej(L\) where the areas inside and outside 
corotation cannot communicate superimposed on the density distribution of the real particles of this Jacobi constant value. 
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Figure 21. Similar to Fig. 20 but for the manifolds of the PLi^ and —1:1 families and for Ej 



-1120000 just above Ej(Li). 



asymptotic curve, and we have integrated all these points 
for 10 iterations (Fig. 22a) 15 iterations (Fig. 22b) and 20 
iterations (Fig. 22c). We observe the formation of character- 
istic rays along the asymptotic curves when integrated for a 
long time. 

In Fig. 23 we see the picture of the same Poincare sur- 
face of section after the integration of 100 initial conditions 
of test particles for 200 iterations in the 2D approximation 
of the system (black points) . These points are concentrated 
close to the unstable asymptotic curves of the families PLi, 
PL2. In Fig. 23 we also plot the curves corresponding to 
zero inertial energy E — (thick black curves). The ana- 
lytic relation z = f(z) (for y = 0) describing these curves is 
found as follows: By subtracting Eqs. [2] and [6] for E — and 
x = x = y = we find 



where 



Hence 



y = ±\/2Ej + Q 2 z 2 - 2V(0, 0, z) - z 2 



-2V(0,0,*) 



site 



(9) 



(10) 



Q p yz + QpZ 2 = 



(8) 



This last equation gives two symmetric curves, one with 
z < and the other with z > 0, shown by thick black curves 
in Fig. 23. Orbits starting below the lower thick black curve 
escape to infinity without intersecting the E = curve. An 
example is given by the blue dots, numbered 1,2,3,..., which 
correspond to the successive iterations of an initial condition 
near the point 1, having E > 0. These iterations are located 
along rays. The corresponding orbit escapes from the galaxy 
following a spiral path, in the rotating configuration space 
(see inserted orbit in Fig. 23). On the other hand, if an 
initial condition has E < 0, then the successive iterations 
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Figure 22. The asymptotic curves of the unstable periodic orbit PL2 in the phase space (z, z) (for y = and y > 0) for the 2D 
approximation of our model and for Ej = —1120000 close to Ej{L\) after (a) 10 iterations (b) 15 iterations and (c) 20 iterations of a 
small initial segment along the unstable directions of the manifold. 



(red dots) lie on rays again, but between the two E = 
curves, and at the same time they form "bell- type" curves 
(Contopoulos and Patsis 2006). Different bell- type curves 
correspond to different values of the energy E. The value 
of E along an orbit remains roughly constant far from the 
galaxy (left part of Fig. 23) forming a bell-type curve, but 
it changes abruptly whenever the orbit comes close to the 
main body of the galaxy (small absolute values of z) and 
then it follows a new bell- type curve. The changes of E are 
irregular, due to the chaotic character of the orbit. After 
several changes of its value, E becomes positive, and then 
the orbit escapes to infinity following a spiral-like path. Note 
that this path corresponds to an open hyperbolic curve in 
an inert ial frame of reference. 



Finally, Fig. 24 presents the fraction of the particles 
starting close to the 3:1, 2:1, PLi,2 periodic orbits (see 
Fig. 13-15) that stay located inside a radius of 2rhm as 
a function of time in Thmct- This figure testifies the phe- 
nomenon of stickiness along the manifolds originating from 
the unstable periodic orbits and gives a measure of the 
chaotic diffusion of the system. We observe two different 
rates of diffusion: (a) the first one lasts for 50 to 100 Thmct 
where about the 20%, 25% and 45% of the orbits around the 
2:1, 3:1 and PL 1,2 periodic orbits respectively have been dif- 
fused outwards and (b) a quite slower diffusion starting after 
& lOOThmct- In general, orbits with initial conditions closer 
to lower order n of resonance n : 1 present a slower diffusion 
rate. During the first diffusion period, the spiral structure is 
viable, because the flow of chaotic material coming from in- 
side corotation stay confined close and along the asymptotic 
manifolds of the unstable periodic orbits and support the 
spiral structure of the galaxy. This time period corresponds 
to about 10 rotations of the bar. After that and during the 
second period of slow diffusion where chaotic orbits move 
outwards, spiral structure is not clearly observed in the sys- 
tem. 
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Figure 23. The phase space 2, z (for y = and y > 0) of 
a 2-D approximation of our model for a Jacobi constant Ej = 
— 1120000 > Ej(Li), where the area inside corotation can com- 
municate with the area outside corotation. The thick black curves 
correspond to zero inertial energy. An initial condition near the 
point 1 and inertial energy E > has successive iterations (blue 
dots) along rays and escape to infinity following a spiral like orbit 
(inserted orbit). On the other hand if we take initial conditions 
having E < the successive iterations (red dots) lie along rays 
and along bell type curves at the same time. 



6 CONCLUSIONS 

In this paper we study the main factors that affect the for- 
mation and the longevity of the stellar spiral arms of barred 
spiral galaxies. We examine the role of the apocentres, of 
the pericentres and of the velocity minima as regards the 
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Figure 24. The percentage of the orbits starting close to the 3:1, 
2:1, PLi ? 2 unstable periodic orbits (see Figs. 13-15) that stay 
located inside R = 2r^ m as a function of time in T hmct . 



morphological features of a galaxy. We emphasize the role 
of the asymptotic orbits, emanating from the unstable peri- 
odic orbits, in forming different segments of the spiral arms 
corresponding to different values of the Jacobi constant in- 
side and outside corotation. We show that there is stickiness 
along these asymptotic orbits which affects the longevity of 
the spiral arms, which finally fade away, but after more than 
10 bar rotations. Below we present and discuss our main re- 
sults. 

(a) An objective way to find the maxima of density is by 
finding the loci of the velocity minima on the rotation plane, 
since the particles spend most of their time there. These loci 
are: 

(i) The apocentres for the particles with Jacobi con- 
stants below El 1 that are allowed to move only inside 
corotation. These apocentres, of regular or chaotic orbits, 
support the shape of the bar. The regular orbits form the 
main body of the bar, while the chaotic orbits form an 
envelope of the bar. 

(ii) The pericentres for the particles with Jacobi con- 
stants below El 1 that are allowed to move only outside 
corotation. These pericentres support the outer parts of 
the spiral structure of the galaxy. 

(iii) An elliptical locus corresponding to the local maxi- 
mum of the effective potential V e f / for the particles with 
Jacobi constant Ej > Ej(L^) that are allowed to move 
both inside and outside corotation. This locus is only 
partly realized in the distribution of the real A/"-body par- 
ticles, because an important percentage of particles with 
Ej > Ej(L±) have minimum velocities v m in that do not 
differ significantly from the mean velocity < v yz > all 
along their orbits. Therefore, this locus is not reflected in 
the global configuration of the particles although it is re- 
flected in the sub-population corresponding to low v m in 
values compared to < v yz >. 



For Jacobi constants Ej{L±) > Ej > Ej(Li) some of the loci 
of pericentres and apocentres that are close to low values of 
the velocity are correlated to the inner parts of the spiral 
structure. 

(b) For each value of the Jacobi constant there is a set of 
stable and unstable periodic orbits. It is well known that the 
domain of phase space associated with the stable periodic 
orbits supports features similar to the morphologies of these 
orbits. 

On the other hand in the present study we emphasize the 
role of the unstable periodic orbits and the manifolds em- 
anating from them. Although the unstable periodic orbits 
themselves do not support the spiral structure, the asymp- 
totic orbits that start near the unstable periodic orbits sup- 
port the spiral structure. As long as there is a particle popu- 
lation near the unstable periodic orbits of the main families 
of orbits (e.g. -1:1, -2:1, PL U PL 2 , PL 4 , PL 5 , 4:1, 3:1 and 
2:1) the spiral structure remains clearly visible. We show 
that the superposition of orbits with initial conditions close 
to the periodic orbits of only three main families (-1:1, 2:1, 
3:1) all along their characteristics are able to reconstruct all 
the main morphological features of the galaxy. 

(c) For Jacobi constants above Ej(Li) chaotic diffusion 
leads to escapes. The diffusion is fast (for all the families) 
for a time interval lasting for about 10 rotations. After- 
wards, the diffusion rate slows down but the spiral structure 
is marginally discernible. For the families inside and close to 
corotation (e.g. 2:1, 3:1, 4:1, PLi, PL2) the diffusion rate is 
smaller for orbits closer to the center. On the other hand, 
chaotic orbits having initial conditions outside corotation 
present stickiness at the resonances —1:1 and —2:1 and 
support the outer parts of the spiral structure. 

The stickiness of the chaotic orbits along the unsta- 
ble asymptotic curves, is ca lled "stickiness in chaos" 
(|Contopoulos &; Harsou la 2008), as it is not due to islands 
of stability or cantori. The shape of these asymptotic curves 
in the phase space ("rays") corresponds to the shape of the 
spiral structure in the configuration space. 

In our model, the phase space outside corotation is not 
bounded. Therefore, chaotic orbits above and outside coro- 
tation will finally acquire positive energy and escape from 
the system, along spiral orbits (in the rotating frame of 
reference). The successive iterations of chaotic orbits with 
negative inert ial energy, in the phase space, lie on bell- type 
curves, each corresponding to the same energy. However, the 
energy changes when a star approaches the main body of the 
galaxy. Once a chaotic orbit acquires positive energy it has 
successive iterations that lie on rays that extend to infinity. 
This escape will happen, in general, after times considerably 
longer than the age of the Universe. 

Finally, we note that our model is an adiabatic approx- 
imation of a real Af-body system. Such a study is capable 
of providing the main trends about the dynamical mech- 
anisms of the formation and the longevity of the various 
structures. However, the dynamical behavior of the real sys- 
tem is definitely more complicated since it shows an evolu- 
tion (decrease) of the pattern speed (although this decrease 
is rather small). Another issue has to do with the evolu- 
tion of the density and consequently of the potential. In our 
"frozen" model the amplitude of the spiral perturbation in 
the density decreases with time. The change of these two 
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parameters affects the phase space structure. This implies a 
different network of unstable periodic orbits and also differ- 
ent asymptotic orbits with different stickiness around them. 
An extended study of the real evolving system requires the 
use of dynamical methods that take into account the evolu- 
tion of both the pattern speed and the potential. We have 
proceeded in this direction and we will present our results 
in a forthcoming paper. 
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